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A recently developed technique to determine the order and strength of phase transitions by extracting the den- 
sity of partition function zeroes (a continuous function) from finite-size systems (a discrete data set) is generalized 
to systems for which (?) some or all of the zeroes occur in degenerate sets and/or (ii) they are not confined to 
a singular line in the complex plane. The technique is demonstrated by application to the case of free Wilson 
fermions. 
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1. INTRODUCTION 

Lattice regularization of a quantum field the- 
ory renders the system a statistical mechanical 
one. The issues of phase transitions and their 
properties become, therefore, of central impor- 
tance. While a true phase transition can only 
occur for a system of infinite extent, the non- 
perturbative computational approach to lattice 
field theory and statistical physics accesses only 
systems of limited size. The partition function of 
such a system can be written as a polynomial in 
an appropriate temperature-like or field-like vari- 
able and the complex zeroes of such a polynomial 
encode all of the information on the behaviour of 
standard thermodynamic quantities. 

Traditional statistical mechanical techniques 
involving partition function zeroes are mainly 
confined to analysing the zeroes closest to the real 
axis, as these are the strongest contributors to the 
critical or pseudocritical behaviour. However, a 
full understanding of the critical properties of the 
infinite-size system requires knowledge of the den- 
sity of zeroes too. It has long been known that 
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the density of zeroes for a finite-size system would 
provide a lucrative source of information but a re- 
liable technique for the extraction of this quantity 
from numerical data proved elusive. 

Recently, however, some of us have succeeded 
in providing just such a technique 0j2|. The gen- 
eral idea is to focus on the integrated density of 
zeroes rather than directly on the density itself, 
which, for a finite system, is a series of delta- 
functions. This new method has been seen to 
be quite reliable and robust and is applicable to 
phase transitions of the temperature-driven 
and field-driven types and to transitions of first 
and higher order. For a comparison between this 
and other approaches, see ||. 

For systems hitherto analysed, the zeroes of 
the partition function had two special properties, 
which seem to be common to the bulk of standard 
models encountered in statistical physics. These 
are (i) the zeroes are all simple zeroes (zeroes of 
order one) and (ii) they lie on a curve called the 
singular line, which impacts on to the real axis 
at the phase transition point. 

The question now arises as to the generality 
of the techniques developed in |l],||]. Here, we 
show how the methods can indeed be extended 
to systems for which some or all of the zeroes 
occur in degenerate sets and/or they are not con- 
fined to a singular line, but instead form a two- 
dimensional pattern in the complex plane. Such 
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two-dimensional patterns of zeroes have been ob- 
served in various lattice field theory and statis- 
tical physics models ||,||. Here, we demonstrate 
the extended technique by application to the case 
of free Wilson fermions. The zeroes of this lattice 
field theory display the two new features we wish 
to address. 

The partition function for a lattice of finite 
extent, L, is Zl(z) oc Yij i z — z j{L)), where z 
is an appropriate coupling parameter. In the 
case where the zeroes, Zj, are on a singular 
line impacting on to the real axis at the crit- 
ical point, z c , they can be parameterised by 
zj — z c + fj exp (if). The density is then de- 
fined as g>i(r) = L~ d J2j <K r ~ r j{L))- The cumu- 
lative distribution function of zeroes is Gl(t) = 
Jo 9L(s)ds, which is j / L d if r 6 (r^rj+i). At a 
zero one assumes the cumulative density is given 
by the average 

G L (r j ) = (2j-l)/2L d . (1) 

In the thermodynamic limit, for a second-order 
transition, the integrated density is, in fact §, 

G=oW«r 2 - , (2) 

where a is the usual critical exponent associated 
with specific heat. Standard finite-size scaling 
emerges quite naturally from this approach JjJ . 

2. GENERAL DISTRIBUTIONS OF ZE- 
ROES 

A departure from smooth linear sets of zeroes 
was found in 1984 when it was shown that for 
anisotropic two-dimensional lattices there can ex- 
ist a two-dimensional distribution (area) of zeroes 
PJ|. Since then, a host of systems have been dis- 
covered with this feature. A common character- 
istic of all two-dimensional distributions of zeroes 
is that they cross the physically relevant real axis 
at only one point, in the thermodynamic limit, 
corresponding to the phase transition. 

For such two-dimensional distributions, the 
density of zeroes in the infinite-volume limit has 
been shown to be 

g(x,y)=y 1 -— /(^) , (3) 



where (x, y) give the location of zeroes with the 
critical point as the origin. Integrating out the 
^-dependence yields 

9(y)= g(x,y)dx ocy 1 " , (4) 

where x\ and X2 mark the extremities of the dis- 
tribution of zeroes at a distance y from the x axis. 
Integrating again, to determine the cumulative 
density of zeroes at the point r in the y-direction, 
yields an expression identical to (Q). Thus, the 
strength of the transition, as measured by a, can 
be determined by similar methods to those pre- 
viously used. Rather than counting the zeroes 
along the singular line, one now counts them up 
to a line y = r within the two-dimensional domain 
they inhabit. 

The second new feature we wish to accommo- 
date is the existence of degeneracies in the set 
of zeroes. If a number of zeroes coincide, Gl, 
as defined in ([j]), is multivalued and is no longer 
a proper function. A more appropriate density 
function is determined as follows. Suppose, in 
general, that Zj — Zj + ± = . . . = Zj +n ^i are rt-fold 
degenerate. It is easy to convince oneself that the 
densities to the left and right of an actual zero are 

r M _ / j + n - 1 for r e (r j+n _i,r J+n ) ( , 
UL[r) ~\ j-1 for re fo-i,^). 1 J 

The density at the n-fold degenerate zero, r,-, is 
again sensibly defined as an average: 

Gr ^i) = ^(i + ?- 1 ) • ( 6 ) 

This is the most general formula for extracting 
the density of any distribution of zeroes and deals 
with two-dimensional spreads and degeneracies. 

3. APPLICATION 

We wish to test the above technique in a sit- 
uation with the two new features of interest - 
namely a two-dimensional set of degenerate ze- 
roes. The free Wilson fermion model is an ideal 
testing ground and its zeroes in two dimensions 
are easily generated exactly § . The zeroes for a 
system of size L — 40 are depicted in Fig. 1 in the 
complex z — 1 /2k plane. Here k is the hopping 
parameter. 
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Figure 1. The partition function zeroes for free 
Wilson fermions on a 40 x 40 lattice, where z = 
1/2k is the inverse hopping parameter. 
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Figure 2. The distribution of zeroes for free 
Wilson fermions. The symbols □, O, x, o and 
o correspond to square lattices of length L = 20, 
50, 100, 150, and 200, respectively. 



We generated zeroes for a wide range of lat- 
tices, and their distributions, as given by (0), are 
plotted in Fig. 2. The data collapses onto a uni- 
versal curve which goes through the origin, indi- 
cating that (0) is indeed a suitable form for the 
density of zeroes and indicating the occurence of 
a phase transition. Fits close to the origin yield 
an exponent compatible with the expected value, 
a = 0. The error estimates appropriate to such 
a fit are non-trivial and we leave their discussion 
to a separate publication . 

4. CONCLUSIONS 

A new method to extract the (continuous) den- 
sity of zeroes from (discrete) finite-size data has 
been extended to deal with the case of two- 
dimensional distributions of zeroes and systems 
in which the zeroes occur in degenerate sets. The 
method has been demonstrated in an application 
to the free Wilson fermion model and seen to be 
capable of direct determination of the strength of 
the phase transition as measured by the critical 
exponent a. 
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